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ABSTRACT 



29 30 

The excitation function for the Si (p>^^)P reaction has been 
experimentally determined in the proton energy range from 1420 kev 
through 2160 kev. Unreported resonances in this reaction have been 
identified at 1975, 2033, 2075, and 2117 kev with estimated errors 
of + 4 kev. Gamma spectra at these resonances are obscured by 
contaminant radiation but suggest complex decay schemes. Double 
angular correlation measurements have been made at the 1470-kev re- 
sonance indicating a resonance level spin of 2+, or less probably, 
2°o The level at 4.50 Mev is shown to probably have a spin of 2+, 
or less probably, 3+ rather than 0+ as has been predicted. 
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lo INTRODUCTION 



As a result of successive refinements and extensions of the collect 

tive model of nuclear structure^ it has been possible to extend its ap« 

plicability into the region of intermediate mass nuclei (20^A^40)» By 

considering the outer nucleons to be strongly coupled to the nuclear core, 

25 25 28 29 

excited states in A1 > Mg , A1 , and Si have been interpreted as col“ 

1-4 

lective rotational excited states with considerable success. More re- 

31 30 

cently states in P and Si have been identified with collective vi- 

5 6 

brational states of a weakly coupled collective model. ^ The experimental 
collection of information on the states of intermediate mass nuclei serves 
to verify the theoretical models, to determine the parameters which enter 
into descriptions of the models, and to permit the extensions and refine- 
ments to theory necessary for a more adequate description of nuclei. Of 
particular relevance to the present investigation is the prediction based 

on the weakly coupled collective model applied to A “ 30 nuclei that the 

30 6 

4.50-Mev level of P has a spin-parity of 0+. In the present work the 

experimental determination of this value has been the object of double 

angular correlation measurements performed on the gamma rays observed at 

29 30 

the 1470-kev resonance in the Si (p, iT)P reaction. 

29 30 

Determinations of resonance energies in the Si (PstJ)? reaction 
have been carried out by numerous investigators since 1941.^ The pro- 
ton energies which have been employed in these investigations extend over 
the range 300 to 1950 kev. 

In addition to detecting resonances and measuring the energies at 
which they occur, the allied problem of analyzing the associated gamma-ray 
spectra has been the object of an increasingly large amount of research 
commensurate with improvements in gamma-ray spectroscopy techniques and 
instrumentation. 

Investigations of the gamma rays which arise from deexcitation of 

30 

the excited states of P were begun in 1954 by Endt, Kluyver, and Van 
1 6 29 3 0 

der Leun. They studied the Si (p,^)P reaction at resonances oc- 
curring at proton energies of 414 kev and 326 kev and assigned spin, 

parity, and isotopic spin to the ground state and the first excited state 
30 

of P premised on gamma- transition probabilities, shell model consider- 
ations and the isotopic spin multiplet foirmalism. In 1956 Broude et al.^^ 
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confirmed the assignments to the first excited level by measuring the 

angular distribution of the gamma^rays at the 414“kev resonance. They 

30 

also detected five new states in P by examining gamma radiation 
arising from the same reaction at resonances occurring at proton ener« 
gies of 696 kev and 737 kev, 

17 18 

Investigations by Lee and Mooring and by Endt and Paris , who 
32 30 

studied the S (d,a)P reaction^ showed the first excited state to be 

in reality a close doublet. The latter investigators were able to accu- 

30 

rately determine the energies of 30 excited levels of P up to excitation 
energies of 5.8 Mev by magnetic analysis of the alpha-particle groups from 
this reaction, 

19 

In 1958 Van der Leun and Endt reinvestigated the four aforemen- 
tioned resonances, studying the decay schemes by analysis of gamma-ray 
spectra, coincidence measurements, and angular distribution measurements. 
From this work came decay schemes, branching ratios, and assignment of 

spin-parities and isotopic spins to the four resonance levels and to five 

30 

other excited levels of P , including the two states of the above men- 
tioned doublet. During 1962 the gamma-ray spectra associated with the 

20 

proton resonances at 1330 kev and the resonances at 1375 kev and 1505 



kev were investigated 



21 



Proposed decay schemes were advanced for these 

30 

resonances and spin-parity assignments made to three levels of P , One 

29 30 

of the more comprehensive investigations of the Si (Ps/*)P reaction 

22 

is that of Baart et al. , who studied six resonances and proposed decay 

schemes for them. Firm spin-parities were assigned to six low-lying 

levels and tentative assignments were made to four more levels. The 

23 

most recent investigation of the reaction is that by Moore, As a re- 
sult of this work decay schemes were proposed for resonances at 1302, 
1470, 1505, 1640, 1686, 1748, and 1772 kev, and spin-parities were de- 
termined for three resonance levels and two intermediate levels. The 

30 

present knowledge of the quantum states of the lower P levels re- 
sulting from the above investigations is summarized in Fig, 1, Two com- 
pilations which summarize important experimental results of the above 

24 

mentioned investigations are the Nuclear Data Sheets and the compi- 
lation of information on the energy levels of nuclei with Z = 11 to 

25 

Z = 20 by Endt and Van der Leun, 
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The purpose of the present work has been to investigate resonances 
29 30 

in the Si (p,?T)P reaction in the proton energy range from 1950 kev 

to 2150 kev which has not previously been reported on, and to investigate 

the gamma-ray spectrum associated with the known resonance occurring at 

1470 kev. In particular^ the angular distributions of gamma rays at the 

latter resonance have been measured in a double angular correlation ex« 

30 

periment to obtain information on the 4.50“Mev level in P . Additionally, 

some information has been obtained from examination of spectra taken at the 

19 

higher resonances although contamination by F hampered this work. The 
following sections of this paper describe the experimental equipment, pro- 
cedure, and data analysis employed, and the results which have been ob“ 
tained from the investigation. 
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II. DESCRIPTION OF EQUIPMENT 



1. Experimental Layout 

Fig. 2 indicates the general arrangement of equipment used in this 
experiment. Fig. 3 is a photograph of the target area showing the cor° 
relation table tilted into a vertical plane. 

2. Accelerator 

The 2“Mev Van de Graaff horizontal electrostatic accelerator of the 
UoS. Naval Postgraduate School, manufactered by the High Voltage Engi- 
neering Corporation, was used in this work as a source of protons. The 

26 

radio frequency ion source is of the type described by Moak. A proton 
spin«resonance fluxmeter for accurately measuring the analyzing magnetos 
field was installed furing the course of the work. An error signal from 
the fluxmeter was used in a feed back loop to the power supply of the 
beam analyzing magnet to hold the field strength automatically at a set 
value. This feature simplified the maintenance of a given resonance 
energy for long bombardment times without "’slipping off*'. 

A current integrator of the vacuum tube electrometer type was used 
to measure beam current. High stability plastic insulated capacitors 
were used to collect the charge. The voltage across the capacitor in 
use was applied to the grid of a 954 acorn tube, whose output fired a 
Schmidt trigger circuit at a preselected voltage level. The trigger 
output actuated other circuitry to discharge the capacitor and advance 
a Veeder-Root cycle counter after each cycle. 

3. Target Assemblies 

Two types of targets and target mounts were developed and used in 
this experiment. These are shown in Fig. 4. For resonance and coinci*- 
dence measurements, which were made in the plane perpendicular to the 
beam, the target backing consisted of a one^half ‘-i n.-diameter flanged 
copper hemisphere held against a flat copper boss on the mount by a 
threaded retaining ring. This assembly was designed with the beam di° 
rection as a symmetry axis so that it can be used for triple correlation 
experiments in the perpendicular plane with no perturbations from target 
geometry. For the double correlations, measured in the horizontal plane 
containing the beam., target backings were 0.046 cm thick copper strips 
clamped in a copper mount. Both types of backings and mounts were vacuum 
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FIGURE 2 . EXPERIMENTAL LAYOUT 




Figure 3 Target Area 
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FIGURE 4. TARGET ASSEMBLY DETAILS 




plated with a layer of gold thick enough to stop incident protons. For 
both geometries the target mount was fitted inside one end of a one-in, 
diameter tubular chambei ^ the axis of which coincided with the beam path. 
Both target assemblies were designed so as to keep mass in the target 
vicinity to a minimum. Competing with this requirement was the need to 
cool the target in order to prevent any appreciable target material eva« 
poration. Maximum possible beam intensities were desired so as to take 
advantage of the better counting statistics which accompany high yields. 
Still another requirement of angular correlation measurements is a small 
target spot to approximate a point source of gamma rays. A beam colli- 
mator consisting of two five-mil tantalum discs separated by 0.4 in. was 
mounted in a brass holder 38 in. from the target. A 0.096“in. diameter 
collimating aperture was followed by a 0.098“in. hole in the second disc 
for the purpose of stopping protons scattered from the edge of the col- 
limating disc. A cylindrical lead shield surrounded the collimator in 
order to attenuate any gamma rays produced here. 

As is described in Chapter Illg target material was evaporated di- 
rectly onto the gold-plated copper backings to permit good heat conduction. 
Initially, a hollow copper target backing through which coolant could be 
circulated was tried, but was abandoned for a simpler solid copper back- 
ing cooled by conduction through a finned copper extension leading to an 
external Dewar flask. Liquid air in the Dewar proved to be convenient 
as a coolant at high beam powers andj, below about ten watts, target power 
dissipation by forced air cooling of the target extension fins proved 
adequate . 

After unsuccessful trials of Luclte and glass target tubes, aluminum 
was finally adopted for this purpose. An auxiliary target assembly with 
a quartz window target was used to accurately center the proton beam spot. 
The target chamber tube, which was insulated from ground at both ends, 
was maintained at a negative potential of 300 volts telati\^e to ground in 
order to suppress secondary electrons from the target. In order to pre- 
vent recoil protons from exciting aluminum resonances, the target tube 
was lined with one-mil gold foil. 

The hemispherical copper target backings for measurements in the 
plane perpendicular to the beam were needed in some quantity, so were 
"mass produced" in a small hand forge made for this purpose. A surface 
which needed only minor polishing and cleaning before gold plating was 
obtained in this manner . 
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A frequently troublesome feature of (p> ^ ) investigations is the 
gradual accumulation of carbon on the target as a result of cracking 
of pump oil vapor by particle bombardment. To reduce this difficulty, 
advantage was taken of an 8“in. glass cross which was mounted between 
the target diffusion pump and the target chamber as shown in Fig. 5. 

A special liquid air cold trap was constructed and installed In this 
cross to serve a dual purpose. First, it functioned as an ordinary 
cold trap with a large surface of copper vanes suspended from the 
coolant reservoir to provide an optically dense path from the dif- 
fusion pump thus permitting the vacuum to be routinely maintained close 
to 10 ^ Torr. The second feature of the trap was a l°in. diameter 
copper pipe piercing the vanes which the proton beam traversed. At 
high vacuum the proton beam functions effectively as a pump, trans- 
porting gases and vapor encountered in its path. When the beam passes 
through the restricted cold pipe for 6 in. before entering the target 
chamber, it was reasoned that its drifting retinue of condensable vapor 
contaminants would be largely trapped by the cold inner wall. The re- 
sults from use of the trap support this idea since only minor discolor- 
ration was observed on targets compared with that which occurred when 
the trap was not used, and no increase in bombardment energy was found 
necessary as targets aged. 

4. Shielding 

To reduce the radiation level at the target resulting from brems- 

strahlung and X-rays caused by electrons accelerated toward the high 

voltage accelerator terminal, a considerable amount of lead shielding 

was erected. A 2-in. “thick lead wall, 5 ft. high and 4 ft. wide, was 

placed between the energy control slits and the beam collimator. In 

addition, two lead screens 1/4-in. thick and each 5 ft. by 4 ft. were 

bolted onto a caster supported steel framework constructed for this 

purpose. The mobility of the screens permitted theii placement to 

minimize the radiation level at the target and also allowed easy re- 

moval for access to the target area. Although the above-described 

shielding offered some reduction in natural background radiation, 

40 

(notably the 1.46-Mev K gamma ray) it was not feasible to screen 
specifically against background without either increasing scattered 
target radiation at the detectors or utilizing an excessive amount of 
lead. 
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Figure 5 Cold Trap Area 




5. Correlation table 

A correlation table shown in Figo 6 was built for this experiment 
to support the gamma-ray detectors and to permit their accurate positioning 
with respect to the target o A steel base plate mounted on casters supported 
a hollow magnesium colum.n upon which the table was mounted. Screw jacks in 
the base plate were extended when the apparatus was in place to hold the 
table in position and to allow accurate vertical adjustment. The table it“ 
self consisted of a movable graduated aluminium disc which was mounted con- 
centrically on a fixed aluminum support plate. Slides mounted on the disc 
and on the support plate held mounts for the two gamma-ray detection as“ 
semblieso The mounts were positioned in the slides by threaded rods which 
allowed the detectors to be set accurately at distances from 1/2 in. to 
6 in. from the target position. The detection assembly m.ounted on the disc 
could be positioned in angle with respect to the fl^ed detector in 12 min. in- 
crements of arc by means of an adjusting gear meshing with gear teeth on 
the periphery of the disc. The entire table assembly was attached to the 
magnesium support column by a hinge so that it could be tilted into either 
a horizontal or vertical plaoe^ permitting measurements to be taken in 
either the horizontal plane containing the proton beam or in the plane per- 
pendicular to the beam. The first named arrangement is designated herein 
as geometry A and the latter as geometry B. 

6. Gamma-ray detectors 

The gamma-ray detection assemblies were supported with their center- 
lines 8 3/4 in. above the table. They each consisted of a 3 in. x 3 in. 
cylindrical sodium iodide (thallium activated) crystal integrally mounted 
with a Dumont 6363 photomultiplier tube supplied by the Harshaw Chemical 
Company. In the earlier phases of this workj, similar 2 in. x 2 in. crystals 
mounted with Dumont 6292 tubes were used. The resolution of the 3-in. cry- 
stals, measured after the signal had traversed a 12-microsecond delay line 

137 

input to the pulse height analyzer, was 97o for the 0.662-Mev Cs gamma ray. 

A 6AK5 cathode follower was constructed into each detector mount to serve 
as a preamplifier. 

7. Electronics 

The electronic equipment and circuitry for processing signals from 
the detectors were planned so as to permit considerable flexibility in ex“ 
perimental applications. A block diagram of the basic layout is shown in 
Fig. 7. 
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Figure 6 Correlation Table 
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FIGURE 7, BLOCK DIAGRAM OF ELECTRONICS 



The system is composed of an RCL 128“channel pulse height analyzer 
to record pulse height spectra^ a slow-fast-sum coincidence analyzer, and 
scalers for recording various outputs. 

It was found that requirements for the detector output signal by 
several components, including a monitor oscilliscope, tended to overload 
the preamplifiers and to permit unacceptable interaction of subsequent 
stages with each other through the coupling of their inputs. To avoid this 
problem a bank of six cathode followers was constructed, A bussing arrange- 
ment made it possible to feed the grids of the followers from either one of 
two busses or from an individual input, so that up to six pieces of equip- 
ment could be driven by the same input pulse with no interaction and with 
practically no signal loss. Each follower was provided with an individual 
output gain control. 

The output of either detector as desired was also fed directly into 
the pulse height analyzer through a simple resistance isolation network. 

From the resistance network the signal passed through a terminated 1350-ohm 
delay line to allow time for formation of gating pulses needed for coinci- 
dence spectra. Several methods of injecting a signal into the analyzer 
were tried inasmuch as this equipment was designed to accept an input only 
from a photomultiplier tube connected directly into the preamplifier as«= 
sembly provided. This was not an acceptable practice because there existed 
no convenient means of signal extraction for other purposes, and because 
the preamplifier utilized a high-voltage supply which required running the 
photomultiplier cathode at a high negative potential with respect to the 
tube shield with the attendant probability of tube damage. The procedure 
adopted as most successful vas to apply the signal from the delay line to 
the anode pin of a lA-^pin tube base which was connected directly into the 
analyzer preamplifier. It was also necessary to incorporate a gate pulse^ 
forming network into the analyzer in order to meet the gating pulse speci- 
fications required by this equipment. 

Gamma-ray spectra were obtained by gating the analyzer with the dis*= 
criminator output pulses from the linear amplifier fed by the same de- 
tector, This procedure excluded a weak precursor pulse which was trans- 
mitted by the delay line inptt, by not gating the analyzer on until a few 
microseconds before arrival of the delayed pulse. No significant loss in 
detector resolution resulted from leaving the delay line in the input for 



15 



ordinary gamma-ray spectra, and it was found quite convenient to be able 
to switch from these singles spectra to coincidence spectra without change 
in analyzer calibration. 

For the observation of coincidence spectra, the analyzer was gated by 

the slow-fast coincidence output. The basic slow-fast coincidence system 

27 

assembled for this experiment has been described elsewhere. In this 

system the output from an Atomic Instrument Company slow coincidence unit 

accepting up to four inputs was combined in a diode mixer with the fast 

coincidence unit output to obtain the slow-fast coincidence output signal. 

28 

The fast unit, based on the circuit of Draper and Fleischer, detects co= 
incidences between amplified pulses from the scintillation detectors with 
a resolution time as low as ten ns. For pulses spanning a large amplitude 
range, however, more dependable operation was obtained when the resolving 
time was increased to about 80 ns. This was due to an observed amplitude 
dependence of the shape of the leading edge of the amplified pulse, to 
which the fast unit is sensitive. 

Several modifications have been made to the coincidence system which 

do not appear in the previous description. The diode mixer unit for com^ 

bining the slow and fast coincidence outputs has been modified by inserting 

a Schmidt trigger between the adding stage and the output stage in order to 

clip the sum pulse and improve output-pulse characteristics. The fast unit 

has been modified by providing a variable delay line in the input in order 

to permit accurate compensation for differences in time delays between the 

detectors and the coincidence unit. Finally, there has been added to the 

system a provision for obtaining sum-coincidence spectra as first described 
29 

by Hoogenboom. In this system the pulses from the two detectors are 

carefully matched in amplitude calibration, then summed in a pulse adding 

circuit. These summed pulses are analyzed in a differential discriminator 

set to accept pulses corresponding to the excitation energy of a resonance 

level. This output is then used as a third input to the slow coincidence 

unit so that its output represents coincidence between two completely ab° 

sorbed gamma rays comprising a double cascade from the resonance level. 

This signal is then used to trigger the pulse height analyzer gate. The 

spectrum obtained from pulse height analysis of one detector output thus 

consists only of full energy peaks due to radiations which are members of 

a double cascade. The pulse height adder which has been constructed for 

this purpose employs the circuit used by Kantele and Fink in a somewhat 

30 

different type of spectrometer. The adder output was analyzed by a con- 
ventional linear amplifier and discriminator. 
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The scalers used to record single channel counts had dead times of one 
microsecond and those recording coincidences had dead times of approximately 
eight microseconds. Separate cathode^f ollower outputs to drive scalers were 
constructed for all coincidence stages so the relatively low scaler input 
impedances would not degrade coincidence pulses » The number of scalers pro“ 
vided, which exceeded the requirements for the present experimental measure^ 
ments, permitted a continuous check on the performance of the coincidence 
system. 

A double pulser with associated pulse shaping networks and an os“ 

cilloscope formed an integral part of the equipment and were used for 

all routine checks and adjustments. An additional slow coincidence unit 

31 

of variable resolving time , which was constructed as part of the system, 
was used to check on the operation of the Atomic Instrument Company unit. 

Control circuit relays were provided so that the pulse height analyzer 
and all scalers could be started and stopped simultaneously by the current 
integrator on“off switch. 
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HI. EXPERIMENTAL PROCEDURE 



1. Target preparation 

The target material used in this experiment was isotopically enriched 

silicon in the dioxide form purchased from the Stable Isotopes Division of 

Oak Ridge National Laboratory. Relative percentages of the three stable 

29 28 

silicon isotopes were 91.85 Si , 7.85 Si , and 0.30 Si . This material 
was evaporated onto prepared target backings in a vacuum evaporation ap= 
paratus under a chamber pressure of 10 to lO” Torr. 

Beam power of as much as 20 watts was delivered to a target spot area 
of about 6 square millimeters during bombardment. Dissipation of this heat 
by 5-mil tantalum target backings, which had been used in previous work at 
this laboratory, was found insufficient to prevent local overheating of the 
target with an accompanying rapid loss of target material. This fact, and 
the fluorine contamination of the surface of all sheet tantalum examined, 
led to employment of the heavier target backings of copper described in 
Chapter II. While vacuum plating these copper backings with 60,000 X or 
more of gold (99.999 percent purity) to completely stop 2-Mev protons, it 
was found that the plating tended to crack and spall off when the copper 
backings were maintained at ambient temperature during evaporation. This 
effect was caused chiefly by tensile stresses which have been found to be 
as much as 850 Kg/cm in 2000 R gold films deposited on copper. This 
trouble was overcome by holding the copper backings at 400°C during eva= 
poration by means of a resistance heater which was constructed for this 
purpose. The heater was made in the form of a 4-in. cylindrical section, 
to which the target backings were clamped. The cylindrical heater section 
was suspended in the vacuum coating chamber so that its axis coincided 
with the evaporant source. Heating the backings in this manner prior to 
evaporating both gold and Si 02 also served to remove any residual volatile 
contaminants leaving a clean surface to receive the coating. 

Evaporating boats were fabricated from 5-mil tungsten ribbon and, when 
cleaned by preheating before use, proved to be satisfactory for the eva- 
poration of both gold and Si02. Temperature of the Si02 evaporant was 
monitored with an optical pyrometer to maintain an indicated evaporating 
temperature of 2130 at this temperature 4.0 mg of Si02 evaporated in 
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about two minutes. The temperature is not corrected for departure of 
the evaporant from an ideal black body inasmuch as reproducibility and 
temperature control were desired rather than an accurate measurement of 
true surface temperature. During the first five to ten seconds of all 
evaporations the evaporant was intercepted by a shutter, since the first 
part of an evaporant to boil off is known to carry a higher percentage 
of any contaminants present. Following evaporation, cooling of the 
heater and targets was accelerated by admitting pure helium to the chamber 
through a liquid air trap. 

In practice, 2 gm or more of gold were evaporated at a distance of 
5 cm from the target backings and exposed parts of target assemblies to 
produce the desired plating thickness. 4.1 mg of SiO^, were evaporated at 
the same distance to give targets approximately 5°kev thick to 2=Mev pro- 
tons for use in the excitation function runs. For the targets used to 
obtain gamma-ray spectra, and for angular correlation runs at 1470 kev, 

5.2 mg of SiO^, were evaporated to obtain targets 8-kev thick to 1470-kev 
protons . 

Considerable care was used in an attempt to reduce the contamination 

19 o 

by F . Evaporation boats were preheated to 2200 C before use, targets 

were stored in a vacuum (until this was found not to be of value), target 

chamber pressure was maintained at 10 ^ Torr during runs, and possible 

sources of contamination in the vacuum systems such as Teflon were a- 

19 

voided. In preliminary bombardments at the F resonance energy of 1373 kev, 

the gold-plated copper target backings, which were washed in either pure 

ethanol or menthanol and rinsed in distilled water just prior to use, were 

19 

found to be freer from F than backings of either Ta or W which had been 

cleaned by various methods. In spite of these precautions all targets 

19 

which were prepared proved to be contaminated to some degree by F 
2. Excitation function measurements 

29 30 

The excitation function of the Si (p>JT)P reaction has been 
measured by recording single gamma^ray yield and gamma-gamma coincidence 
yield at accelerator energy settings over a range of proton energies from 
1420 kev to 2170 kev. The 2 in. x 2 in. Nal(Tl) scintillation detectors 
were used for this purpose. The output of one of the detectors was counted 
to obtain the gamma-ray yield at each proton energy, and the number of 
events in which gamma rays were detected simultaneously in each crystal was 
recorded by the use of the slow-fast coincidence system adjusted to a re- 
solving time of 25 ns. 
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The two detectors were positioned in geometry B at angles of 0° and 

180^ relative to the vertical direction and at a distance of 0.6 in. from 

the beam spot. Integral discriminators in both channels were set so that 

counts would be recorded for gamma rays with energies exceeding 0.65 Mev, 

a setting chosen to exclude the 0.51-Mev annihilation radiation arising 

30 

from the positron decay of P , while accepting the full energy peaks of 

the 0.68 and 0.71 Mev radiations which deexcite the first two excited 
30 

states in P . The 7-microampere proton beam was collected by the cur- 
rent integrator which closed a beam shutter and stopped counters after 
collection of 84 microcoulombs at each energy step of approximately 3 kev. 
Background counts were made every 50 kev with the beam intercepted by a 
shutter in the beam tube outside of the shielded target area. 

An auxiliary experiment was conducted to determine the energy cali- 
bration for the excitation curves. In this experiment the frequency of 

the NMR signal was measured to within one part in 10 ^ as the 992-kev 

33 

aluminum resonance was traversed using a freshly prepared target about 

10 kev in thickness. Runs with mass 2 beams were also made to check this 

calibration at a generator energy of 1984 kev. Immediately after cali- 
29 

bration, the Si resonance energies above 1970 kev and the energy of 
the resonance at 1506 kev were carefully measured in a similar manner. 
Using these points as references it was then possible to calibrate the 
entire excitation curve in terms of magnet current, which was determined 
to a precision of one part in 10 ^ at each energy step. This procedure 
gave a better calibration than that obtained using the fluxmeter logging 
dial readings, since calibration curves for converting these readings to 
magnetic field could not be read with adequate precision. 

3. Measurement of saturation activity 

29 30 

In order to definitely identify resonances due to the Si (p,y)P 

reaction, measurements were made at the resonances above 1970 kev to 

determine the half life of decay radiation following bombardment. In 

29 

these measurements Si targets were bombarded for ten minutes, then re° 

moved to a small cave of 2-in. lead brick where counting was started 

within 60 sec after cessation of bombardment. The positron annihilation 

30 

radiation following beta decay of P was detected by a 3 in. x 3 in. 
crystal. Total counts above the noise threshold were recorded by the 
pulse height analyzer operated in a gross counting mode. In this mode 
the channel in which all counts are recorded is automatically shifted 
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ahead after a preset counting time. After counting for 3.1 sec per 
channel for a 6.46-minute total counting time, the target was removed 
from the cave and a background run made with the analyzer in the sub- 
tract mode. 

Difficulty which was experienced in obtaining gamma^ray spectra 
29 

recognizably due to Si at the upper resonances as a result of masking* 

radiation from the F^^(p,o: 8^)0^^ reaction caused some concern that the 

peaks observed in the excitation curve might have been due to other con“ 

30 

taminants. Observance of the 2,5-min half life of P following bombard- 

ment could not entirely exclude this possibility since it is conceivable 
30 

for P to be formed in a nonresonant direct capture process o For this 

reason it was desirable to compare the activity of the target observed 

after bombardment at energies corresponding to the peaks in the yield 

curves to the activity observed after bombardment at off-peak energies. 

In an experiment to do this, a proton bombardment time of ten minutes to 

obtain 94 percent of saturation activity was followed by a ten-minute 

counting period during which the pulse height spectrum was obtained. 

The target was then allowed to decay for an additional ten minutes before 

bombardment was resumed at the next energy. The number of counts observed 

in the 0*51-Mev annihilation radiation peak then served as a measure of 
30 

the amount of P formed at each energy. 

4. Gamma-ray spectra 

All gamma-ray spectra were obtained using 3 in. x 3 in» Nal(Tl) 
crystals with no special shielding nor collimationo For most spectra 
the pulse height analyzer usable range extended from about channel 7 
to channel 127 corresponding to gamma-ray energies from 0.7 to 7.7 Mev« 

For observation of annihilation radiation a higher gain was used so that 
these channels corresponded to energies from 0.18 to 1.92 Mev. The 
analysis of spectra for intensities was performed with a digital computer 
as described in Appendix 2. The particular method finally used to sub- 
tract background radiation from spectra varied according to the resonance 
being considered. When possible, backgrounds taken slightly above or be- 
low resonances were subtracted from spectra. This was not very successful 

19 

for the resonances above 1950 kev where F contamination gave rise to the 

F^^ (p,a 8 )0^^ reaction. This reaction is strongly resonant at 2030 kev 

34 

with a width at half maximum of 120 kev which effectively blankets the 
higher proton energy range considered in the present work. The 6-7 Mev 
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gamma rays occurring in this reaction are present in relative proportions 
which vary considerably with bombarding energy so that of f ^resonance back° 
ground subtraction is not accurate. Bombardment of blank targets which 
had been processed identically to SiO^ targets but using an empty eva- 
poration boat gave a background measurement which was also rather unre« 
liable because of individual variations among targets. A method of sub° 
tracting the background measured with the accelerator off was also used, 
in which case background radiation of specific energies due to contaminants 
was subtracted out in the spectrum analysis. 

In order to circumvent radiation due to target contaminants efforts 
were exerted to obtain siiri-coincidence spectra in the energy region above 
1950 kev as a means of getting at decay schemes. The sum-coincidence 
spectra were, in general, measured with the slow-fast coincidence require- 
ment also imposed, and with a coincidence resolving time of 110 ns. These 
runs were usually made in geometry B with counters positioned 1.6 in. from 
the beam spot at a 90° angle relative to each other. A 5/16-in. lead sheet 
between the counters served to minimize coincidences due to scattering or 
pair production in one crystal with the lost radiation being captured in the 
other crystal. The differential discriminator following the pulse adder was 
set with its lower level at the resonance energy and was adjusted to pass 
pulses in a band about 0.6 Mev in width. Settings were made by adjusting 
the discriminator so that the observed count rate of accepted pulses was 
equal to half the repetition rate of a test pulse. The proper test pulse 
amplitude was obtained by setting the pulser gain so that counts fell into 
the pulse height analyzer channel corresponding to the desired energy. 

Settings of integral discriminators in the two singles channels were adjusted 

to 0.62 Mev in a similiar manner. The pulse height analyzer was calibrated 
1. 3 7 60 

with Cs and Co sources and with the 5.31-Mev gamma ray which dominates 
29 30 

the Si (p, JT)P spectrum at the 414-kev resonance. It was found that the 

Dumont 6363 photomultiplier tubes experienced marked gain shifts as a function 

of count rate amounting to several percent, an effect which has been attri- 

35 

buted to changes in the dynode surfaces . Attempts to allow this shift 
to stabilize before commencing counts were only partially effective. The 
overall result of this phenomenon was to vitiate the results of a large 
number of sum-coincidence spectra. 
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5. Measurements at 1470 kev 

30 

Since the 4«50“Mev level in P has been shown to be most strongly 
excited at the 1470“kev resonance, this proton energy was a logical choice 
for angular correlation measurements which could lead to determination of 
the spin*=parity for the level. The resonance is rather weak so it was de- 
cided to take five sets of data to improve the counting statistics and to 
permit a check of the consistency of results. Each set of data consisted 
of five spectra taken with the movable detector set to angles of 0°, 30°, 
45°, 60^, and 90° relative to the proton beam. In geometry A with the 
3 in. X 3 in. movable detector positioned with its face 4.125 in. from 
the target, runs were made using a 7 -microampere beam for counting times 
of approximately one hour at each angular setting. The same target was 
used for each set of runs and the order of the angular settings was varied 
in different sets so as to minimize differences in the gamma-ray spectra 
due to any contaminant buildup or target erosion. The fixed 3 in. x 3 in. 
detector was positioned at 270^ relative to the beam and 4.125 in. from 
the target for use as a monitor counter. Each run was terminated when 
30,000 monitor counts had been recorded. A 2 in. x 2 in. detector was 
also used as a secondary monitor counter. The integral discriminator in 
each monitor channel was adjusted to a valley in the gamma^ray spectrum 
at 4.75 Mev so that the effect of gain shifts in these channels would be 
minimized. Before and after each run the discriminator levels were checked 
and adjusted with the aid of a Co^^ source which was precisely and repro- 
ducibly located relative to the counters. For this check the gain in the 
monitor channels was increased by a factor of 4 by means of the gain mul- 
tiplier switch so that the discriminator setting was equivalent to 1.14 Mev. 
At this setting, near the 1.17“Mev peak in Co^^, the count rate observed 
with the Co^^ source was quite sensitive to the discriminator setting so 
that accurate adjustment was possible. 

The unfolding of gamma^ray spectra and the determination of angular 
distribution coefficients was programed and analyzed on a CDC“1604 computer 
as described in Appendix 2. The radiations which have been used as spectro° 
meter calibration curves in this analysis are listed in Table 1. These 
spectra were measured in experimental geometry A, backgrounds were sub- 
tracted, and any remaining contaminating gamma rays were subtracted off. 

Each of these monoenergetic spectrograms was plotted, hand smoothed, and 
redivided into 525 channels to reduce any interpolation errors in the 
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TABLE 1 

Gamma Rays Used for Spectrometer Calibration 



Ej(Mev) 


Radiation Source 


E (kev) 
P 


0.662 


Cs'^' 




1.63 


90 

Na (p,Ct ^)Ne 


1252 


2.37 


12 13 


460 


2.75 


24 

Na 


“ 


3,51 




1698 


4.43 




675 


5.31 


c-29, ;^.„30 

Si (p, if)P 


414 


6.135 


„19. j,..^16 

F (p,a (jr)o 


1374 


8.06 


^13, y.,„14 

C (p, )N 


554 
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subsequent analysis. Energy calibration of each spectrum was made by 
equating distance in channels between the pair peak and the photopeak 
to 1.022 Mev. For the three lowest energy calibration curves, for which 
pair peaks are non-existent or poorly defined, an equivalent calibration 
was computed based on the pulse height analyzer calibration as determined 
in an auxiliary experiment. 

All spectra were read out of the pulse height analyzer memory photo- 
graphically. These data were then transcribed to punched cards and printed 
listings of the cards were carefully proofread against the photographs. 

Gain drifts did not prove troublesome at the lower count rates ob“ 
served during correlation runs at 1470 kev. To compensate for differences 
in calibration of the pulse height analyzer between successive runs, which 
were of the order of one percent, each spectrum was self “calibrated by the 
peaks at 7.00 Mev and at 1,45 Mev, 
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IV EXPERIMENTAL RESULTS 



1. Excitation Function 

29 30 

The excitation function for Si (p, if )P as measured by gamma-ray 
yield is shown in Fig. 8, and as measured by gamma^gamma coincidence 
yield in Fig. 9. The energies at which resonances are observed are listed 
in Table 2, together with some of the results which have been obtained at 
proton energies below 1900 kev by other investigators. 

In arriving at the uncertainty of 0.207o kev for the resonance energies, 
consideration has been given to errors arising from instability and inhomo- 
geneity of the field of the analyzing magnet, the proton beam spread, pos° 

sible surface contaminants, calibration drift, and the 0.5-kev estimated 

33 

error for the 992-kev aluminum calibration point. 

Particular attention was given to the resonances above 1900 kev which 

have not previously been reported. The results of determinations of the 

half life of the post-bombardment activity of the targets at these energies 

are given in Table 3. The measurements are seen to agree with values for 

30 

the half life of P reported in the literature, which range from 2.51 + 

24 29 

.02 min to 2.62 + .02 min, so that radiative capture by Si is clearly 

demonstrated. Table 4 shows the intensity of positron annihilation 

radiation observed after saturation at selected bombardment energies. 

These intensity values show distinct maxima at 1975, 2033, 2075, and 2117 

30 

kev and show that P is formed in a resonant reaction rather than by a 

non-resonant process. The peaks in the gamma-ray yield curves at these 

energies can therefore be positively ascribed to resonant proton reaction 
29 

with Si 

The yield ofQoSll-Mev radiation after bombardment at energies between 

resonances is probably due to a combination of off^resonant reaction with 

29 17 

Si and to the 66°sec positron activity of F following the non-resonant 

X )F^^ reaction. This hypothesis is supported by half life measure^ 
ments of post^bombardment activity of the targets at two inter-resonance 
points at 1955 and 2007 kev« At those energies *'half -lives*' of the order 
of 1.6 + .5 min are observed with the count rate 2 min after bombardment 
at 1955 kev a factor of ten lower than the count rate 2 min after bombard- 
ment at 1975 kev. 
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FIGURE 8. GAMMA-RAY YIELD 
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TABLE 2 

29 30 

Resonances Observed in the Si (p, (T)P Reaction 

E from 1450 to 2170. Values in kev. 

P 



This work 


Green 


12 

Seagondollar 


10 

Tsytko 


Wild^^ 


Remarks 


1470 


1479 


1479 




1470 




1506 


1515 


1512 




1505 




1637 


1648 


1643 


1635 


1638 


Si^s 


(1646) 

1667 


1671 


(1653) 


(1647) 

1663 


(1650) 

(1665 

(1670 


1688 


1692 


1693 


1680 


1685 




(1726) 










weak 


1752 


1752 


1752 




1749 




1775 


1777 


1775 




1772 




(1803) 


1811 






1808 




1860 


1857 






1852 





(1933) 

1975 

(2007) 

2033 

2075 

2117 

(2159) 



weak 



Si^° 



not identified 



+ 0 . 20 % 



+2 



+2 



+3 



Errors 
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TABLE 3 



Half Life of Post “bombardment Decay of 



P 



30 



Bombardment Energy (kev) 
1975 
2033 
2075 
2117 



Half life (min) 
2,42 + .09 
2.45 + .18 
2.51 + .06 
2.47 + .15 
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TABLE 4 



Intensity of Post=boitibardment Annihilation Radiation 



Bombarding 
Energy (kev) 


ic 

Excitation Relative 

Energy in P (Mev) Intensity 


1506 


7.04 19.3 


1950 


1.0 


1975 


7.49 4.1 


1990 


1.0 


2033 


7.55 7.9 


2060 


1.8 


2075 


7.59 12.4 


2100 


1.3 


2117 


7.63 8.1 


2140 


2.3 



Based on Q = 5.582, reference 24. 



31 



2. Gairana^ray Spectra at 1975, 2033, 2075 and 2117 kev Resonances 

Gainma^ray spectra were obtained at the four resonances at proton 

energies of 1975, 2033, 2075 and 2117 kev. Masking radiation from the 

19 

resonant reaction with the contaminant F dominates the upper energy 

region of these spectra so that a detailed analysis of the relatively 

29 30 

weaker radiation from the Si (p, )P reaction was not possible. 
However, the gamma-ray energies and intensities which could be deter- 
mined with reasonable certainty from stripping and least squares analysis 

19 

for intensities are given in Table 5. The 6=7 Mev gamma rays due to F , 
which are omitted in Table 5, were accounted for by inclusion in the 
least squares intensity analysis, since other background subtraction 
techniques did not work well at these energies as has been discussed 
in Chapter III. Also included in Table 5 are energies of sum-coinci- 
dence spectral peaks which were close to energies observed in the 
singles spectra. The good spectral resolution and discrimination 
against background radiation which are obtainable with the sum-coinci- 
dence method proved to be of only nominal assistance in analyzing the 
gamma-ray spectra. This was due in part to the relative weakness of 
radiative decay of these resonance levels in comparison with inelastic 
scattering, and in part due to photomultiplier gain drifts which blurred 
peaks in the sum-coincidence spectra. However, it is apparent from the 
large number of low intensity peaks observed in the sum-coincidence 
spectra at the higher-energy resonances that decay schemes from these 

levels are considerably more complex than decay from the 1752-kev re- 
23 

sonance level, where sum-coincidence spectra were recorded for the 
purpose of comparison. 

The strong 1.28-Mev gamma ray observed at 1975, 2033, and 2075 kev 

30 

is interpreted as deexcitation of the P compound nucleus by proton 

29 25 

reemission to the first excited state of Si at 1.277 Mev. The high 
intensity of this gamma ray is convincing evidence that this inelastic 
scattering process predominates over radiative capture and that the 
128-Mev radiation is primarily responsible for the upper-energy peaks 
in the gamma-ray yield curve. Rather curiously, however, this gamma 
ray does not appear in the spectrum at the 2117-kev resonance. 

28 v' 29 

The resonance at 2075 kev lies near the Si (p,5)P resonance 

25 

reported at 2090 kev which has a half width of 16 kev. The gamma“ray 
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TABLE 5 





Gamma-ray Intensities at Resonance 
Energies of 1975, 2033, 2075, and 2117 


kev 


Ep (kev) 


Singles Spectra 
E ^ (Mev) 


Sum^Coincidence 
E ^ (Mev) 


Intensity (quanta^/*- C) 


1975 


4.52 


4.58 


5 




4.19 


4.07 


2 




3.48 

3.33 


3.40 


2 

2 




3.00 


2.99 


3 




2.44 


2.50 


3 




1.28 


- 


117 


2033 


5.10 


- 


8.6 




1.28 


- 


187.3 


2075 


7.59 


- 


14 




4.70 


4.85 


60 




2.98 


2.95 


55 




2.26 


2.23 


55 




1.36 


- 


303 




1.28 


- 


1360 


2117 


5.10 


- 


2 




2.55 




3 



Error in singles spectra E^ estimated to be 17o 
Error in sum°coincidence Ey estimated to be 2% 
Error in intensities estimated to be 30% 
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TABLE 6 



Relative Gamma-ray Intensities at the 1470-kev Resonance 



E y (Mev) 


Transition 


"k 

Observed 

Intensity 


Decay Scheme 
Intensity 


Ref 23 


7.00 


R-^ 0 


25.4 


25.4 


39.6 


6.29 


R-^ 2 


9.6 


9.6 


6.3 


4.50 


19-^ 0 


14.6 


14.3 


13.2 


4.14 


13-^ 0 


10.7 


10.7 


11.1 


3.98 


R-»-9 


10.5 


10.8 


7.3 


3.45 


not accounted for 


2.0 




4.0 


3.17 


R-^ll 


3.0 


1.5 


2.0 


3.17 


11-^ 1 




1.5 


2.0 


3.07 


R-^12 


19.1 


17.6 


11.7 


3.02 


9-^0 




1.2 


1.6 


2.86 


R^13 


13.7 


13.7 


13.0 


2.69 


13^3 


2.8 


2.8 


- 


2.50 


R-^19 


38.8 


21.6 


20.0 


2.50 


12-^3 




17.6 


11.7 


2.34 


9 1 


17.2 


16.9 


13.4 


1.45 


19'*- 9 


not analyzed 


^ 7.3 


6.8 


1.45 


3-»“0 


II 


20.3 


11.7 


.71 


2->0 


II 


9.6 


21.8 


.68 


l->0 


II 


18.4 





An intensity of approximately 34 is indicated for the 1.45 Mev radiation 
by observations at 90° with .6*' crystal face to target distance. 

Relative intensities shown are also approximate absolute intensities 
in quanta per microcoulomb. 
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Double Angular Correlation Coefficients at 1470-kev Resonance 
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spectrum at this resonance is characterized by a 4o76“Mev radiation 

29 36 

to the ground state of P . It seems probable, therefore, that 

the 4,70“Mev gamma ray seen in the spectrum at 2075 kev is due to the 
28 

Si reaction o The observed intensity is in agreement with this ex- 
planation. Although the presence of 2o94«Mev and 2.24-Mev gamma rays 

30 

in this spectrum suggest a cascade through level 8 at 2.94-Mev in P , 

22 

the lo45-Mev gamma ray known to occur in this cascade is not seen 

and sum-coincidence work has not been able to confirm this possibility. 

29 

It is possible, however, that the gamma rays due to Si reactions could 

be responsible for the anisotropy of the radiation from the 4.76-Mev 

29 37 

level in P which has been observed by Newton, since anisotropy is 

not expected to result from the spin assignment to this level of 1/2+ 

38 

based on proton scattering experiments. 

The remaining weaker radiations seen at the four higher resonances 

29 

are most probably due to radiative capture by Si , but the rather large 

errors, due chiefly to masking contaminant radiation, do not permit more 

than tentative postulation of decay modes. Gamma rays observed at the 

1975”kev resonance suggest cascades through level 19 at 4.50 Mev and 

through level 14 at 4.18 Mev. No conclusions can be drawn for the 

2033-kev resonance. At the 2075-kev resonance the 7.59-Mev radiation 

is readily recognized as a ground state transition but the other gamma 

rays are not identified. Finally, at the 2117“kev resonance, the two 

radiations observed nicely fit a decay through level 5 at 2.54 Mev which 

22 

is known to decay 957c. to the ground state. 

3. 1470-kev Resonance Decay and Intensities 

The gamma-ray spectrum observed at the 1470°kev resonance with a 
crystal face to target distance of 0.6 in. and at an observation angle 
of 90° is shown in Fig. 10. This spectrum has been least squares fitted 
by gamma rays with energies required by the decay scheme proposed in ref- 
erence 23. Reanalyses with alternate energy choices in each instance gave 
a poorer fit in reassurance of the correctness of the decay modes proposed. 
This decay scheme, with intensities as found in the present analysis, is 
shown in Fig. 11. (Heavy dots indicate cascades confirmed in reference 23 
by coincidence spectroscopy.) The branching ratio shown for the 4.50-Mev 
level is that given in reference 22. The assumed correctness of this 
ratio then requires branching ratios of 93:7 for level 9 and 80:20 for 
level 13 in order to fit the intensities of gamma rays which have been 
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Mev 

FIGURE 10. GAMMA-RAY SPECTRUM AT 1470 
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observed. These ratios are in reasonably good agreement with ratios of 

22 23 

89:11 reported for level 9 and 88:12 reported for level 13 

Intensities of gamma rays observed in the angular correlation runs 

are shown in Table 6 along with intensities observed in reference 23 at 

55°. The intensities found in the correlation runs have been averaged 

over 4jt steradians by weighting the mean values for each observation angle 

by the sine of that angle. Because the 1.46-Mev background radiation due 
40 

to K was found to be quite strong relative to the 1.45-Mev 
gamma ray observed in the correlation runs, no attempt was made to analyze 
the gamma radiation of this or lower energy. Although the background 
radiation subtracts out, the statistics of the remaining counts do not 
allow extraction of significant information. The observed intensities 
are seen to agree reasonably well with those of reference 23 except for 
the 7.00“Mev radiation which is found to be considerably weaker in the 
present work. Since the experimental intensities are based upon reported 
values of intrinsic peak efficiencies^^ rather than an absolute spectro- 
meter calibration, this discrepancy is not considered significant. 

4. Angular Correlations at the 1470-kev Resonance 

A. General 

The correlation function coefficients which have been found in this 
experiment are given in Table 7. The values given are weighted means of 
the values found in each of the five sets of data taken. The chi-square 
values for each correlation coefficient have been calculated and show the 
different sets to be consistent, under the acceptance requirement that the 
cumulative chi-square probability be equal to or less than 0.9. Combining 
intensities from each data set, then analyzing for the coefficients gives 
values practically identical to the mean values quoted. This agreement 
to some degree serves as a check that computational errors have not crept 
in. Each experimentally observed spectrum has also been plotted together 
with the corresponding fitted spectrum, then visually examined to determine 
that a reasonable fit was obtained in each case. Fig. 12 shows a typical 
example of the fit obtained for the correlation runs. 

The intensity analysis was performed both subtracting background 
radiation analytically and subtracting of f ^resonance background, with 
essentially the same results being obtained in each case. The standard 
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FIGURE 12. GAMMA RAY SPECTRUM FOR 
CORRELATION RUN 3-5 AT 90® 



deviations in the mean coefficient values have been calculated from 



1 /^ 

i=l ^ 



and have been compared with the values given by 



1 

n(n-l) 



n 



i=l 



(X, - X) 



2 



which agree closely „ 



The most obvious feature of the data is the anisotropy of the 

4.50“Mev gamma ray which is incompatible with a spin of zero for this 

level. It is also noteworthy that the 4p50-Mev gamma ray is the only 

radiation to show a statistically significant P^(cos 0) term (notation 

is defined in Appendix 1). Fig. 13 shows the correlation coefficient 

values for this gamma ray obtained from each set of data. In this 

figure M denotes mean values and C denotes values obtained by com= 

bining intensities from each data set before computing coefficients. 

Fig. 14 shows the fit to the Legendre series both with and without 

terms. Although the existence of the P^ term is not incontrovertible 

it permits a markedly better fit to the data and gives a smaller value 

of the parameter € which serves as a guide to the appropriate number of 

41 

terms which should appear in the Legendre series. 

In order to deduce spin^parity values it is necessary to examine 
the possible particle and radiation mixing ratios for this reaction. 

The notation used herein to refer to these mixing ratios is as follows: 

Proton°partial“wave mixing parameter 

The multipole mixing ratio of an intermediate 
unobserved gamma ray 

The multipole mixing ratio of an observed 
gamma ray 

Channel“spin mixing parameter 

Appendix 1 defines these parameters more fully and describes how they 

enter into the expression for the theoretical correlation function. For 
29 30 29 30 

the Si (p,Y)P and Si (p^xIT)? (intermediate unobserved) reactions, 



Si- 

= 
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RUN I2 345MC I2345MC 



RUN I 2 3 4 5 M C 




HIGHEST TERM 6 ** 1.30 

FIGURE 13. DISTRIBUTION COEFFICIENTS FOR THE 
4.50 mEV gamma ray 
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FIGURE 14. FITS OF COMBINED EXPERIMENTAL DATA TO 
LEGENDRE SERIES - 4.50-MEV GAMMA RAY 
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the correlation coefficients A^/A^ and A,/A^ are describable as functions 

^0 4 0 

of a maximuni of three mixing ratios. Level curves, or contours, of con- 
stant A^/Aq values may be plotted against any two of the mixing ratios. 
Actually the arctangents of the mixing ratios are used as variables for 
convenience since they remain finite. If a third mixing ratio appears, 
changing its value changes these contour plots only to the extent of 
changing the A^/Aq values for all contours by a multiplicative constant. 
Contour plots for this analysis have been programmed and plotted as 
described in Appendix 2. 

In the following discussion it will be assumed that only El radiation 

or E2/M1 mixtures are likely to be observed for transitions from levels 

of several Mev. This simplifying assumption is based upon transition 

42 

probabilities estimated from the independent particle model and is com- 
monly made in analysis of gamma radiation following proton capture in 
light and intermediate nuclei. Although suppression or inhibition of di- 
pole radiation is possible because of the isotopic spin selection rule, 

A 43 

= + 1 for self conjugate nuclei, there are an adequate number of 

30 

lower levels of both T = 0 and T = 1 in P so that the resonance level 
can decay by either El or E2/M1 radiation without violating this rule. 

In analyzing the experimental correlation coefficient values, tables 
from reference 23 have been quite helpful. This reference tabulates the 
maximum and minimum values which may be assumed by A^/Aq and A^/Aq for a 
variety of postulated spin sequences, 

B, Resonance level spin^parity 

The observed A^/Aq terms of about -0,40 for the 7,00, 6,29, and 

3.98 Mev gamma rays, which all proceed to known 1+ levels, are found to 

be most useful in determining the spin-parity of the resonance level. 

Among possible values through 4 for this level, it is found that only 

1+, 2-h, and 2« can lead to the A^/Aq coefficient measured for the 7.00“Mev 

ground state transition. The same conclusion is supported by the A^/Aq 

values observed for the 6,29 and 3,98 Mev radiations. 

If one inspects the contour plot, Fig, 15, for the 1+ — > 1+ transition 

it is seen that a d-s proton-angular -momentum mixing ratio of over 0,90 

is required to produce the observed A^/Aq value for the 7,00-Mev gamma 

ray. Since the relative penetrability of d waves to s waves is about 
-2 

3 X 10 , such a mixture is not a likely choice. The A^/A^ term observed 
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for the 4«50-Mev radiation supports this conculsion since A^^/Aq terms 

do not result from resonance spins smaller than 2. Finally, the value 

of A^/Aq observed for the 2,34“Mev radiation will be shown to be income 

patible with a 1+ resonance level spin. This radiation proceeds to the 

known 0+ level 1, at 0,68 Mev, from the known 1+ level 9^ at 3.02 Mev. 

Level 9 is fed more than 507o from the resonance level and, secondarily, 

by the intermediate 1.48-Mev radiation from level 19. The A^/Aq value 

observed for this primary 3.98-Mev feeder restricts the ranges of 

and ^2 the R— ^9— 5^1 transition. When these ranges are delineated 

on a contour plot for a l+-^l+-^ 0 + transition, it is found that A 2 /AQ 

for the second radiation is restricted to a range between - 0.01 and 

48 

+0.24. By an extension of the formalism for intermediate unobserved 
cascades given in Appendix 1 it is possible to calculate the angular 
distribution coefficients for the 9 1 radiation which result from the 

R->*19-^9“>1 cascade. It is found that, for a 1+ resonance level, 
the minimum possible value of A2/AQ occurs for an assumed spin of 2 for 
level 19. This minimum value has been calculated to be -0.175. Thus 
a 1 + resonance level leads to a minimum possible value for A^/A^ of 
about -0.093, Comparing this with the observed value of -0.24 + ,03 for 
the 2. 34 -Mev gamma ray, 1+ can be ruled out as a possible spin-parity 
of the resonance level. 

Having narrowed the possible resonance spin values to 2+ or 2-, 
the contour plots may be examined for the 2 + 1 + and 2 “->-l+ transitions, 
applicable to the 7.00-Mev radiation to the ground state. In the first 
case, channel-spin mixing may occur while in the second case proton- 
partial-wave mixing is permitted. For the 2+-^l+ case Fig« 16a shows 
the contour plot with the experimental value of A^/A^ for the 7,00-Mev 
gamma ray delineated thereon by cross-hatching. Fig, 16b is a similar 
plot for A^/Aq, The darker shading in Fig, 16b indicates the superposi- 
tion of cross-hatching on both plots and defines the ranges and combi- 
nations of mixing parameters which can give rise to the observed coeffi- 
cient values. For S 3 (i.e. Ml radiation), is constrained to be 

less than ‘^tan 50^. This must not be taken too seriously, however, 
since the A^/Aq contours are widely spaced at the top and bottom boundaries 
of the permitted region, so that a very slight increase in the error for 
the experimental A^/A^ value would allow S ^ to range up to 90^, The small 
region at '^tan(-67^) and ^^^tan 49^ also fits the observed corre- 
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lation. The plots for the 2 case, shown in Figs, 17a and 17b, are 

quite similar to the 2+-^l+ case and show that there is no restriction 
on f-p proton-wave mixing for (El radiation). The two small per- 
mitted regions at tan 67^ may be ruled out since they require an 

M2/E1 mixture. It is seen that an unlikely but small M2/E1 mixture with 
^^^0o026 or greater is required for this transition for the case of 
resonance formation by pure p waves. As will be shown, the 2- resonance 
level assignment is also less probable than 2+ based on mixing ratios 
which are required by experimentally observed coefficients for the 4,50- 
Mev radiation. 

As a result of the above considerations, a spin-parity of 2+ is pro- 
posed for the resonance level with 2- as a less likely possibility. 

C. R— ^13“>0 cascade 

There are only two cascades which allow cross checking of mixing 
parameters as a method of eliminating possible spin sequences. These 
are the cascade through level 19 and that through level 13 both of which 
levels have undetermined spins. The latter cascade will be considered 

first. The existence of an A, /A^ term which has been observed in the dis- 

^ ° 23 

tribution of the 4.14-Mev radiation from level 13 at the 1686-kev resonance 

is evidence that level 13 must have a spin greater than or equal to 2. Of 
the likely possible spins for this level, 2+, 3+, and 4-, only 2+ and 2- 
can result in the negative value of ot>served for the 4.14-Mev gamma 

ray. An examination of the contour plots for the four sequences 2+-^2+->l+ 
for the R ->13-^0 cascade shows compatibility in all cases with the ob- 
served coefficients for the 2,86-Mev gamma ray and the 4,14-Mev gamma ray. 

For the sequence 2+“^2“-^l+ the parameter is constrained to be larger 

than tan 70^ but the other possible sequences are compatible with a full 
range of values for either £ ^ or S ^ as the case may be. The only finding 
is that, for all cases, the R->13 transition must be relatively pure di- 
pole radiation. The observed coefficient values for the weak 2.69-Mev gamma 
ray, which also is believed to deexcite level 13, are also compatible with 
all four possible spin sequences with reasonable choices for S 

D. R->19*->0 cascade 

Since no requirements on parameters applicable to this cascade are 
imposed by transitions already considered it becomes necessary to consider 
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the R->'19->0 sequence with no prior restrictions on the deltas involved 
We find, furthermore, that no information is forthcoming from the R— >19, 
2„50-Mev gamma^ray correlation coefficient since this radiation is indis“ 
tinguishable from the 12->3, 2o48-Mev gamma ray. 

A survey of maxima and minima of coefficients shows, if we 

allow 2- as a possible resonance level spin, that the following sequences 
for the R-^19-->0 transition are consistent with the experimental co- 
efficients for the 4.50-Mev radiation: 



2 + 1 + 1 + 
2 - 1 + 1 + 
2 + 2 + 1 + 
2 - 2 + 1 + 
2+ 3+ 1+ 



2 + 1 - 1 + 
2 - 1 - 1 + 
2 - 2 - 1 + 
2- 3- 1+ 
2+ 3- 1+ 



It has been found in the analysis that a better fit (see Fig. 14) 

to the data is given by the series including a termo Then, provided 

this term is not spurious, which is not to be excluded as a possibility, 

the list can be reduced by eliminating sequences not giving rise to 

A, /A^ terms o 
4 0 

The following possibilities then remain: 

2 + 2 + 1 + 

2 - 2 + 1 + 

2+ 3+ 1+ 

2+ 3- 1+ 

2- 3- 1+ 

The last two transitions shown require E3/M2 radiation so, adhering to 
the hypothesized El and E2/M1 radiative transitions as practical candi- 
dates, the spin of the 4»50-Mev level is limited to either 2+ or 3+. Of 

the three remaining possibilities, the 2 >2+->l+ sequence is perhaps 

least likely. Figs, 19a and 19b show that this sequence for R->19“>0 
requires an f to p wave mixing parameter of over 0,10 with a large E2/M1 
mixing ratio of nearly 12,0, For a 2+->3+->l+ sequence. Figs. 20a and 
20b show that an E2/M1 ratio of over 1,7 is required for the 4,50-Mev 



47 



radiation in order to fit the data. Because large E2/M1 mixing para= 
meters are not commonly seen in the lighter nuclei, this sequence, though 
more likely than 2 ~ 2 + ->-l+, does not appear as probable as a 2 +-^ 2 +-^l+ 
sequence, as will be seen. For the 2+-^2+-^l+ sequence Figs. 18a and 
18b show that for an assumed channel“spin mixing ratio of zero, an E2/M1 
mixing ratio of Oo27 will fit the observed coefficients for the 4.50-Mev 
gamma ray with a considerable latitude permitted for ^ Alternate 
choices for will also permit acceptable choices of mixing parameters 
for this sequence though in general $2 and will be required to be 
larger . 

Summarizing, level 19 is found to have a spin of 2+ or 3+ under the 
assumptions which have been made, with 2 + being more probable. 



E. Comparison with other results 

Concurrently with the present work, an investigation of triple 

angular correlations at the 1470-*kev resonance has been undertaken at 

40 

Wright-Patterson Aeronautical Laboratories. The angular distribution 

of the 7.00-Mev gamma ray was also examined in this experiment. An ^ 2 ^^ 

value of -0.38+<»02 was observed for this gamma ray in good agreement with 

the value shown in Table 7. Preliminary results of this experiment have 

indicated restriction of the resonance level spin to a value of 1 or 2 

and show that a zero spin for the 4,50~Mev level is not compatible with 

6 2 

such a restriction. The preliminary analysis of the triple“Correlation 
experiment is seen to be in agreement with the present work. It is to be 
hoped that the complete analysis of this work which is now in progress, 
will lead to a firmer indication of the spin-parity of the 4.50°-Mev level 
since the 2+ (3+) value proposed herein depends upon the somewhat debatable 
A 4 /A 0 observed and upon the hypothesis that only El or E2/M1 radiation 

is present. 



F. T = 1 character of level 19 

The prediction of 0+ for the spin-parity of level 19 is based on the 

30 

premise that this level is the T = 1 analog of the 3.80-Mev level of Si 

42 

in the isotopic spin multiplet formalism. A spin of 0 has been suggested 

30 

for the 3o80-'Mev level in Si based on observations of cascade gamma rays 
30 30 63 

following the Si (PsP’)Si reaction. A spin of 0+ results also from 
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calculations made by Thankappen and Pandya of collective vibrational states 

in a weakly-coupled collective model of the nucleus.^ The spin-parity of 
30 

level 19 in P is thus expected both on experimental and theoretical grounds a 
The finding of non 0+ spin-parity for this level then appears to be good evi- 
dence that either the level is complex (e.g. a doublet) or that it is not the 

30 

isotopic spin multiplet analog of the 4,80-Mev state in Si as has been ac= 

30 

cepteda Some evidence has been reported that the 3,80-Mev level in Si is 

a doublet, based upon observations of protons from the reactions 

28 30 30 30 59 44 

Si (t,p)Si and Si (p>P*)Si . ^ It is possible, therefore, that 

30 

the analogous 4.50-Mev level in P is a doublet, although no other evi“ 

dence so suggests a The other possible explanation is that a nearby level 
30 

in P , such as that at 4.42 Mev, may be the T = 1 level rather that the 

4.50 Mev level. This does not appear unlikely, since the T = 1 character 

attributed to the 4.50-Mev level is a probable rather than a firm assignment 

32 30 

which was deduced from the investigation of the S (d,a)P reaction by Endt 

18 30 

and Paris. It is based upon a correction of the 3.80-Mev level in Si 

for the Coulomb energy and neutron-proton mass difference which yields a pre- 

30 

dieted energy of 4.47 Mev for the excitation energy in the P nuclide. The 

assignment is strengthened by the fact that the alpha-particle group from 

32 30 30 

S (d,a)P corresponding to a residual excitation in P of 4.501+ .010 Mev 

has only about half the intensity of neighboring groups. This has been inter^ 

preted as being due to operation of the selection rule conserving isotopic 

spin in this reaction. 

30 

On the basis of the energy predicted for the T = 1 level in P the 

4.42-Mev level is seen to be nearly as likely a candidate as the 4.50-Mev 

level, particularly since the energy prediction is not expected to be exact. 

Similarly in the T conservation argument, the observed lower intensity of 

the alpha group for the 4.50-Mev level can hardly be considered conclusive 

evidence of the T = 1 character of the level, since some T = 0 levels are 

associated with alpha-particle groups of even lower intensity (e.g. the 

30 

3.84-Mev level). Unfortunately, the 4.42-Mev level in P is not known to 
be excited in radiative capture, so that spin-parity of this level will pro- 
bable have to be determined from a different reaction. 

5. Summary of Results 

Excitation curves for the Si^^(p,K)P^^ reaction have been obtained over 
the proton energy region from 1420 through 2160 kev. Four previously unreported 
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resonances have been detected and confirmed at proton energies of 1975, 

2033, 2075, and 2117 kev with experimental uncertainties of +4 kev. 

Gamma°ray spectra at these resonances show that inelastic scattering pre= 

dominates over radiative capture, except at 2117 kev where no inelastic 

scattering is evident » Some information concerning the radiative decay 

schemes at the 1975, 2075, and 2117 kev resonances was obtained by spectral 

analysis with the assistance of sum-coincidence spectrao 

The gamma-ray spectrum at the 1470-kev resonance has been analyzed 

23 

and found to be in agreement with the decay scheme proposed by Moore. 

The angular distributions of gamma rays observed at this resonance have 
been measured in a series of five experimental observations, each con- 
sisting of gamma-ray spectra obtained at five angles relative to the pro° 

ton beam. Analysis of this data, under the assumption that El and E2/M1 

30 

radiation deexcite P energy levels of several Mev, has shown that the 
1470-kev resonance level has a spin°parity of 2+ or, less probably, 2-. 

The spin°parity of level 13 at 4.14 Mev has also been found to be limited 
to 2+ or 2-. Level 19 at 4.50 Mev was found not to have a spin=parity of 
0+ which had been predicted. Under the above-mentioned assumption, an 
apparent A^/k^ term observed in the angular distribution of this radiation 
leads to the conclusion that the spin-parity for level 19 is 2+ or, less 
probably, 3+. 
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FIGURE 16a 
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2 + I + A 4 /A 0 
FIGURE 16b 
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2-I + A4/A0 
FIGURE 17b 
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CONTOUR MULTIPLIER = .75-l-.25sin(2tan'S4 +90°) 

FIGURE 18a 
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tan'Si 

2+2+I + A4/A0 

CASE SHOWN FOR 84= 0 
CONTOUR MULTIPLIER= -.127 +.635sin 

(2tan'84-90®) 



FIGURE 18b 
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2“2 + I'**A4 /Ao 
FIGURE 19b 
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2+3+1+ A 4 /Ao 
FIGURE 20a 



60 




61 



V ACKNOWLEDGEMENTS 



I wish to thank my advisors. Professor HoEo Handler and Professor 
E,A. Milne for their encouragement, guidance, and advice offered during 
the course of this worko I also am most grateful to Dr, L.W» Seagondollar 
and Dr, R.A. Moore of the University of Kansas for pointing out the signi“ 
ficance of the 4o50-Mev level spin determination and to Dr, Moore for fur- 
nishing results of his experimental work. Among the many other people who 
have helped make this experiment possible, I wish to thank Mr. R.C. Moeller, 
Mr, M.A. Brillhart, and Mr» R.A, Garcia for their indispensible contributions 
to the experimental equipment, 

I also wish to acknowledge the help and understanding of my wife Lois 
during the years consumed by this undertaking, I wish, finally, to express 
appreciation of the cooperation and confidence of the Civil Engineer 
Corps of the U.S. Navy which have made possible my absence from more normal 
duties in order to complete this work. 



62 



VI APPMDICES 



Appendix 1 - Theoretical Description of Double Angular Correlations 
A . General 

There has been abundant literature coverage of the theory of angular 
correlations of successive nuclear radiations since Hamilton first pre- 
dicted the effect and developed theoretical ejq^ressions for successive 
emission of quanta. A comprehensive review of the subject has been 
published by Biedenham and Rose and recent articles useful in exper- 
imental applications have been published by Gove and by Biedenham. 

The following notation is essentially that of Biedenham. In all cases a 
primed symbol represents a larger value than the same symbol without a 
prime. 

In the case of a (p,^ ) reaction in which the direction of a primary 
gamma ray deexciting a resonance level is observed relative to the direc- 
tion of an incident proton beam^ the angular correlation is referred to 
as a double correlation. The function describing the relative probability 
of the gamma decay as a function of the angle between the beam direction 
and the emitted gamma ray can be written as an expression in the even 
Legendre polynomials s 



W(e) - 2 A (1)A (2)P (cos©) 

TV - A " 



0 



, even 



( 1 ) 



where 
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and F^(LL'j^J) : (-1)^^“*^“^ [(2L + l)(2LVl)(2J+l)j 

(LL*1 -llvO)W(JJLL*j 2 / j^) 

The index h is related to the index i by: 

h s i^ - i + 1 



(Ll'i -1\Z/0) and W(JJLL* j j^) are respectively the 

Clebsch-Gordon coefficients and the Racah coefficients which are commonly 
used to describe systems involving coupling of angular momenta. The F 
functions, for which an extensive tabulation is available, are defined 
in terms of a Clebsch-Gordon and a Racah coefficient by Biedenham and 
Rose. Arguments of the coefficients are angular momenta of the se- 
quence : 



^1 

l^,l; 

J 

Si 






/I’l 




II] 


r 


j 




[hj 










( 2 ) 



= channel spin “ proton spin© target spin 
= interfering proton orbital angular momenta 
= spin of resonance level 
= interfering gamma ray multipolarities 
= spin of final state (not necessarily a ground state) 

= mixing parameter representing the ratio of the reduced 
matrix element for the transition involving l| proton 
waves to that for waves 

= mixing parameter representing the ratio of the amplitude 

f 

of gamma rays of multipolarity L 2 to those of multipo- 
larity L 2 
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( 3 ) 




cos(fc, -L,) . 

L(L+ 1 ) + L (L*+ 1 ) - 2.(X/ + 1 ) 



j^fare the phases at the nuclear 
surface of the particle partial wave functions of multipolarity L and l' . 
The coherent interference between partial waves is expressed by the 

electrostatic potential for charged particles and the phase shift occur- 
ring in "hard sphere" scattering. The 212/^2^2^ terms for photons as 
opposed to particles are equal to unity. In the case where mixing of 
channel spins and occur, W(0) may be expressed as an incoherent 

mixture ; 



) term. The factor ^ includes the coulomb phase shift due to 



is a particle parameter in which ^ and ^ 



W( 9 ) = W^j^(6)+ 52 

c ^ 

where is a mixing parameter expressing the ratio of probability of 

channel spin entering the reaction to that of channel spin 

The maximum order of Legendre polynomials appearing in ( 1 ), ZJ » is 

max 

determined by conditions imposed on the W (abcdj ef) coefficients by con- 
servation of angular momentum which may be expressed by the so-called 
triad relationships. These require each of the triads of values (acf), 
(bdf), (abe), and (cde) to form a possible triangle. Restriction to the 
case of intermediate states of well defined parity (i.e. no interf erring 
intermediate states of opposite parities) leads to the further simplifi- 
cation that may take on only even values. 

In the case of the reaction (p,x^), in which x represents an unob- 
served intermediate gamma ray, the correlation is expressed as an incoher- 
ent mixture of multipoles occurring in the unobserved transition, i.e. in 
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the sequence: 

W(0) cen be expressed: W(0) I 

where: 




(5) 

( 6 ) 



Jq - channel spin 

- interfering proton orbital angular momenta 
= spin of resonance level 

Lj^,L^ = interfering gamma ray multipolarities of unobserved 
radiation 

~ spin of intermediate level 

^ 2*^2 “ interfering gamma ray multipolarities of observed radiation 



= spin of final state 

^2 ” ratio of intensity of the ;uiobserved gamma ray of higher 
multipolarity to the intensity of the lower multipolarity 
gamma ray. 

The correlation function now contains an additional part expressing 
the connection between the two intermediate levels by the unobserved 



radiation as follows: 



'max 



W(e) 2 A^(0)C^(1)A^(2 )P^(cos9) 



(7) 



where the Aj^(i) are the same as in (1) and 



C^(l) = ( -if [(2Ji1-l)(2J2-Hlj] Vji V2? ‘^1'^2^ 



It is seen that for the general double correlation case four mixing 
parameters may occur (barring the further complication of mixed inter- 
mediate states) . 
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In the present investigation the Si29 target nucleus has a ground- 
state spin of 1 / 2 +, with the simplifying result that channel-spin 
mixing and proton-partial-wave mixing may not occur simultanously., The 
number of (S 's is thus limited to a maximum of three. Explicitly^ for 
resonance levels of 14- ,2-, 3+> 4-, etc. orbital angular momentum mixing 
of two partial waves may occur but channel-spin mixing is prohibited. 

For the resonance level spins of opposite parities from those listed 
only channel-spin mixing is permitted. 

When three mixing parameters occur in the general case of inter- 
mediate unobserved radiation (p,x with a spin sequence 




the normalized correlation function can be expressed as 



w(e) = i/AoZap(cos0) 



(9) 



where each A /A^ is a product of the form: 






( 10 ) 



with 



yv - 1,2,3 or 2,3,4 



and 













(It is to be noted that the coefficient A ^ appearing Ln (9) is a conven- 
tional representation of the more complete forms shown in (1) cind (7) . 

A distinction must also be made, of course, between the correlation func- 
tion W(0) and the Racah coefficient W(abcd; ef). 



67 



The coefficients a, b, and c can to expressed as follov/ss 



( V = 

protcn 

wave / s 

mixing 



VI 



2Lo(Lq-^i) p*,(Lo^oJi) 

2Lq(Lo+ 1) -2>(w-H) 




4 [lo(Lq 

Lo(Lo-»- 1 )+ LqW ■*' + 



Fz.UoI^J 0 ^l) ( 11 ) 



2lo(Ln 1) 

2 Lo(L^+ 1 ) » M*'+l) 




unobso 

gamma 

mixing 



/ 

°«/2 = 






^x /3 



gaimna 

mixing 










(-l)'*^[(2Ji+ 1)(2J2+ 1)]^ WCJiX'I^Jgl J 1 J 2 ) 
0 

(“1)*^[(2J2_+ 1)(2J2 ■+■ 1)]^ W(J2^i'L2^J2l'J]L'J2^ 



^l/(^ 2 ^ 2 ^ 2 ^ 

2Fj^(L2L2J2^2^ 



F^(L^J 2 «^ 2 ) 



a 



chan. 

spin 

mixing 






^b /4 



\ ®x 4 



2 Lo(lo + 1 ) - ^ 1 ) 

0 



2 Lo(Lo+ 1 ) 1 ) 






tl - 



- "f 3 - ~f 4 - ^ 
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#4 ^ 



For the (p, J ) case the only changes required are: *^2 “ 

c ^2 “ 0 “ ^2 “ 

The angular correlation function, as an alternative, may be ex- 
pressed in tenns of the coefficients Z and , instead of the F coef- 
ficients* 

f M = (-1)^® z(ViVi> Joi') 



proton 



wave 

mixing 



Ci^l 



2(-l)^0 ZCLqJj^LqJi, jo^) 
(-1)^0 ZCL^^JlL^Ji, JqV ) 




unobs. 

gamma 



- a^2 
< K2 



mixing 






(-1)”^^ W(JiJ2J3^J2l Lj_ Z') 
0 

(-l)’^i W(JiJ2JiJ2J 1^2^) 



Si o 

•*^3 



gamma 

mixing 









®*^3 



(-1)^2 qcqwa* h'"'> 

(-1)«52 Z3^(L^J2L2J2p ^2^^ 



chan. 

spin 

mixing 



» I 

XJk 



1^4 



*^4 



(-1)^0 Z(LoJlVl* 

0 

(“1) ^ Z(LqJ2^LqJj^| JqV) 
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For (p, ir) reactions - Jp = J, 2 ® » 0^ c^2 “ 

* ( The two formalis.'tis (11) and (12) are to some extent incocipatiblec 

. CQ Cl 

Not only may the coherent irixing paranietar phases differ by XI ^ 

but also for = 0 the parameter S| doesn’t represent the same matrix 

eleirent ratio in both notations. The notation above^ (12) j has been used 

for most calculations in the present work. These expressions are in 

agreement with Sharp except that the sign of b^^ ^ (Pj3C^) 

case only, has been changed to agree with the part of Ruby’s correction 

applying to the observed gamma ray,^*^^^^ For the ’’unnatural" radiation 

mixtures (M1/E2, S3 /^i etc.) the sigi of b^ should be reversed. 

This sign convention pennits the use of Moore’s tabulated values 

where the sane convention is used.) 

The Z and coefficients and Racah and Clebsch-Gordon coefficients 

have been tabulated by Sharp et al, Using this tabulation Moore 

has prepared a convenient table of the a, b, and c coefficients applica- 

20 

ble to radiative proton capture by Si 

In order to fit experimental data on angular correlation coeffi- 
cients to theory, it is desirable to first eliminate as many spin se- 
quences as possible on the basis of any evidence which may be available. 

in 

In cases where only two mixing parameters enter, contour plots ^ show- 
ing contours of constant theoretical values of the coefficients 
can be plotted against the inverse tangents of the two mixing parameters 
for all possible spin sequences. If three mixing parameters enter, one 
of them may be held constant for the purpose of plotting the correlation 
coefficient contours against the inverse tangents of the other two param- 
eters. It is seen from (10) that the contour values for different 
values of the fixed parameter differ only by a multiplicative constant. 
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4 






« * - 









V/hen experimental values of correlation coefficients are delineated on 
the contour plot, some spin sequences can be eliminated and the possible 
combinations of mixing parameters are restricted for other sequences. 

The observed values of several correlation coefficients in one decay 
scheme may then fix or highly restrict values of the spin-p 2 irities of 
resonance and intermediate levels, and determine the values of the mix“ 
ing parameters. 

B, Application to 1470-kev resonance 

For the resonance occurring at a laboratory proton energy of 1470 
kev, the phase factor cos' k-fL ,) has been calculated using the 
expression where represents the phase shift due to the 

point coulomb potential andf^ is the "hard sphere” phase shift, 

is defined by exp(2i<T ) s "P (L 1 i<>t ) (I 3 ) 

^ ^ ra + i-i'T^) 



where 



"I 



ZZ'e^ 

“fiv 



which leads tos 

^L + 2-S’l. = 




+ tan“^ 




(U) 



The value for (p ^ defined by 



exp(2i^,) - ° 

where F and G are the regular and irregular solutions to the radial wave 



equation for the coulomb potential. The value for A can be found from 
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the graphs of Sharp et al. (these graphs actually list-^) and em= 
1/0 

ploying R - 1.5A '•'^f, as an estimate of the nuclear radius. These 
graphs have also been used to calculate the relative pentration factors 



Vj/v. 



L-“2 “**■ 



' l<-2 

irCL 

which are useful as a guide to the expected mixing parameter of inter- 

42 

fering partial proton waves. Results of these calculations appear in 
Tables 8 and 9. 



+ G 






( 16 ) 
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TABLE 8 



Relative Penetration Factors for l470“Kev Protons 
d/s f/p g/d 

3.1 X lO"^ 5.0 X lO"^ 1.3 X lO"^ 



TABLE 9 

Coulomb Phase Factors for 1470=Kev Protons 
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Appendix 2 - Computer Programs and Analysis 



Ao General 

In order to analyze experimental data in the form of gamma-ray 
spectra recorded by the 128-channel analyzer , two computer programs 
were developed. The first and more extensive program entitled program 
SPECFIT computes experimental correlation coef f icients^ while program 
CONTOUR is designed to utilize a California Computer Products Model 
565 graph plotter as an output medium to automatically plot contours 
of constant angular correlation coefficient values against the inverse 
tangents of two mixing parameters » Both programs are written in fortran 
language and were run on a CDC-1604 digital computer. 

Bo Program SPECFIT 

The output pulses of a photomultiplier tube viewing a scintillation 
crystal upon which monoenergetic photons impinge do not all have the same 
voltage. Instead, the three fundamental processes of photoelectric ab- 
sorption, Compton scattering and pair production, which occur in the 
crystal with probabilities dependent upon gamma-ray energy, lead to a 
response pulse-spectrum which is characteristic of each particular gamma- 
ray energy. Experimental geometry, crystal size, shielding, various line- 
broadening effects, and characteristics of individual crystal-photomulti- 
plier tube combinations influence the response also, so that response 
functions are in general unique for each spectrometer in a given experi- 
mental arrangement. 

An experimentally observed spectrum is, of course, usually composed 
of a number of, gamma rays, the spectrum being the sum of the constituent 
radiations each multiplied by its intensity. The problem of deciphering 
such an observed spectrum to obtain the energies and intensities of com- 
ponent gamma rays is of fundamental importance in gamma ray spectroscopy 
and has not yet been satisfactorily solved in complete generality. 

Using notation, which is appropriate for the discrete pulse height 
intervals into which multi-channel analyzers sort pulses, and also ap- 
propriate for the discrete number of gamma ray energies composing re- 
sonance deexcitation spectra, one can writes 



m 
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where 



S. = observed number of counts in channel i 
1 



R. . = probability that a photon of energy E. 
originating at a particular source ^ 



will produce a pulse falling in pulse 
height channel i (usually referred to 
as the response matrix) 



I. = the number of photons of energy E. 
^ originating at the source ^ 



m = the number of gamma ray energies which 

occur in the spectrum under consideration 



A rather fundamental and obvious approach to this problem is to 
assume no knowledge of gamma-ray energies is available, divide the gamma 
ray energy range into a number of discrete intervals equal to the number 
of channels, then solve for: 



In fact, the first attempt to unfold experimental spectra in the 

54 

present work utilized a modification of a program which employed an 

55 

iterative technique developed by Schofield to solve the above equation. 
This attempt failed for two reasons which are rather instructive. Firsts 
the analytical expression used to approximate the response matrix did not 
adequately represent the true spectrometer response, and second, the 
technique proved inadequate with experimental counting statistics attain^ 
able to handle the analysis of approximately 100 channels. 

In order to avoid the first problem resort was had to the technique 
of interpolating between experimentally-determined monoenergetic cali- 
bration curves in order to obtain the response matrix. In order to over- 
come the second objection it was decided to use information on the gamma 
ray energies available from preliminary "stripping" of spectra and from 
other sources. Having thus reduced the number of gamma rays to be con- 
sidered, a linear combination of the selected response matrix elements 
was then fitted to the observed spectrum by a least squares process. 

This procedure was followed by another least squares analysis to obtain 




J 



S 



k 
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correlation function coefficients from the gamma ray intensities cal- 
culated for a series of runs at different experimental angles » The 
steps in the computation procedure in program SPECFIT as finally de- 
veloped are as follows: 

lo Load data, (See section D of this Appendix) 

2o Subtract any constant background from calibration spectra. This 
background represents the Compton “tails** of higher energy gamma rays 
contaminating certain of the calibration gamma rays, 

3o Normalize each calibration spectrum to a photopeak height of one 
unit . 

4o Include a straight line in the response matrix if this option is 
selected. This step is essentially a refinement of step 2. 

5o Correct the input and background spectra for dead time of the pulse 
height analyzer. For this purpose the dead time given by the manufacturer 
is 38 + 0.5i microseconds where i is the pulse height channel number. 

This leads to an expression for corrected counts: 



N*. = 
J 



1 - 1/T 



N. 



127 

2 N^(38 + 0.5i)*10 



where N. = the number of counts observed in channel i 
1 

T = counting time in seconds 

6, Calculate energies in Mev corresponding to each channel number of 
the input spectrum, 

7, Adjust energy scale and number of channels of background spectrum 
to match the input spectrum. Multiply by ratio of run times to get the 
proper background to be subtracted from the input spectrum and perform 
the subtraction, 

8, Using the monoenergetic calibration spectra, calculate a matrix of 
response curves for the gamma rays which the spectrum is to be decomposed 
into. This calculation is performed in a subroutine entitled RESPONS as 
follows : 

a. Find the two calibration curves bracketing a gamma ray energy 
desired as a response matrix element. 
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bo Shift the lower-energy calibration curve found in a» so that 
its photopeak falls in the same channel as the higher energy 
photopeak. Fill in bottom of shifted curve as determined by 
the slope between the bottom two unshifted channels. Match 
energy calibration of the lower-energy curve to the upper, 

(For this step, and elsewhere where interpolation is necessary, 
cubic interpolation is performed using the four nearest points 
to the desired value,) 

c. Interpolate linearly between each pair of points of the 

bracketing calibration curves to obtain a response curve for 

the desired gamma-ray energy. (A "curve" in this description 

is meant to imply the vector consisting of a 1 X n array of 

values, n being the number of channels used in the description), 

do Using the input spectrum channel energies calculated in 6,, 

interpolate to get a response curve with channels and channel 

energies corresponding to those of the input data, 

e. Normalize the response curve for the given gamma-ray energy so 

that the ordinate corresponding to each channel represents the 

probability that one photon leaving the source will produce a 

count in that channel. For this calculation the experimental ly^ 

56 

determined intrinsic peak efficiency data of Lazer have been 
used to derive the following approximate expression for intrinsic 
peak efficiency as a function of source distance: 



where A = ,165 + ,0067h 

X =:= ,9699 - ,00918h 

h = source to crystal distance in centimeters 
E = gamma ray energy in Mev 

Assuming an accurately gaussian photopeak, and that the crystal 
resolution is approximately inversely proportional to the square 
root of gamma ray energy, (an assumption which was experimentally 
verified) the ratio of photopeak height to area can be expressed: 



€ = ae'^ 



a 






1 



_ 2 ^ 

WR 




1 

.866WR f~E 
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I 



where R is the usual resolution definition of peak width at 

137 

one half maximum amplitude divided by energy for the Cs 
gamma ray, and W is the spectrometer calibration in pulse 
height channels per Mev. When CX is multiplied by the appro- 
priate solid angle fraction and attenuation coefficient, the 
value obtained represents the probability of recording one 
count in the photopeak channel when one photon originates at 
the source 0 



f o Return to main routine » 

9o Fit the experimental spectrum with a linear combination of the selected 
monoenergetic response curves, the coefficients representing absolute in^ 
tensities because of the method of normalizing the response matrix^ The 
weighted sum of the squares of the differnces between observed channel 
counts and fitted counts is minimized in order to arrive at this fit« 

The normal equations resulting from the minimization are: 



where S. = Observed counts in channel 1 
J 

R., = Response matrix element giving 
^ probability of a count in channel 
j due to one photon of energy 
originating at the source 

w. = Weight assigned to the squared 
^ error in channel j 

I, = Number of photons of energy E. 
which leave the source during 
the counting period 




J 



ij 



j summations are over the channels for 
which a fit is desired 



i summation is over the number of gamma 
ray energies E. used to fit the Ob'- 
served spectrum 



For calculation the are expressed: 





and 
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A standard subroutine (GAUSS3) is used to invert A» Mention should be 
made of the weights w. which have been usedo Normally ^ for a Poisson 

J .2 

distribution the weight is taken as 1/S where the number of observed 

2 

counts in a channel are used as an estimate of the variance S «» 

Ferguson, who has pointed out that this procedure overweights statist!- 
cally low counts, identifies the Poisson distribution variance, which 
rightly is the mean number of counts, with the fitted number of counts 
for each channel » This procedure is adopted here with a modification to 
account for statistical error arising from the subtraction of background 
counts. We have Sl = S. ° b., where S. are observed channel counts cor° 
rected for dead time and b^ are background counts to be subtracted » Then, 

=<$^sp + 

^ S. + (T/T, )b. = s: + (1 + T/T, )b. 

1 ^ b' 1 1 ^ b'^ 1 

where T = counting time for data and T^ = counting time for background 

run. Identifying with the computed number of counts, S V , then gives 

sl* + (1 + T/T,)b. for the channel variances. To avoid overweighting of 
1 b 1 2 

very small values of fitted counts an arbitrary minimum value of S is 

chosen to be two per cent of the average value over the fitted 

channels. Since no fitted count values are initially available, an ini- 
tial weight of one is assigned to each channel and the solution is then 

iterated to convergence, stopping when 

2 

is smaller than 10 

(This series of iterations, usually five in the present work, accounts 
for most of the program running time and can be abbreviated for prelimi- 
nary analyses; see reference 57 in this regard,) is the fractional 

th ^ ^ 

change in the computed intensity of the i gamma ray between successive 
iterations , 

10. The computed intensities are examined for algebraic sign. Since 
only positive intensities have a physical basis any negative intensities 
are discarded and the problem resolved to this point, 

11, For each gamma ray calculate the fitted total number of counts in the 
channels above the channel representing 40 per cent of the photopeak energy. 
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(This is an arbitrary calculation used only for checking results obtained 

in manual unfolding) , The fractional error assigned to these counts is 

set equal to the fractional error in the corresponding intensityo The 

statistical variance estimators of the intensities have been taken to be 

A 1 

the diagonal elements of the inverse of the matrix A in 9o In addition 
to statistical intensity errors there is an error due to gain and discrimi** 
nator drifts in the monitor counts channel which must be accounted for as 
well as error due to eccentricity of the beam spot on the target « These 
fractional errors have been estimated to be 3.3 per cent and 1.2 per cent 
respectively in the present experiment. The three errors are combined 
according to ; 

- 2 VIS’ ^ . 

j ^ 

Two other possible sources of error should be noted here. First, the 
errors in the experimentally determined response matrix, which can con^- 
siderably complicate the error calculations for the intensities, are kept 
to a minimum by obtaining good counting statistics in the calibration 
spectra. The other source of error is gain drift which tends to smear 
spectral peaks. Fortunately, the main effect of such drift seems to be 
an increase in the mean squared errors in the curve fitting process and 
intensity determinations are apparently not too sensitive to the effect. 

12. The intensities and their errors calculated in 11 are read out and 
stored for use in the subsequent determination of correlation coefficients 
and mean values. 

13. The absolute values of differences between observed and fitted counts 
are summed over fitted channels and divided by total input spectrum counts 
over the same range to give a measure of the fit which has been achieved 
in the analysis. This measure is then read out. 

14. The input and fitted spectra are plotted in a graph output. 

15. The above process is repeated for each run in the data set until the 
data for all the experimental angles has been processed. 

16. When the preceding analysis is completed for a set of data a sub° 
routine entitled LEASTSQ may be called, if appropriate, to compute the 
correlation coefficients. For each gamma ray in the spectrum a least 
squares fit of the counts computed at the experimental angles is made to 
the usual series: 

1 

W(6) = 2 A.P.(cos 9), i even 

1 X 
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where are the even Legendre polynomials » The fit is first made to 

Aq + ^2^2 ^0 ^ ^2^2“ (Provision has been made for fitting 

to the series terminating with by simply changing the initial value 

of NP to 4; however, see reference 61 concerning accuracy of P. terms) o 

41 ^ 

The computation, which follows the procedure given by Rose , proceeds 
as follows: 



a» 



b„ 



C o 



do 



e . 



f o 



Table gamma ray absorption coefficients in copper and cor^ 
responding energies in Mev. Correct intensities and sigmas 
for attenuation in copper target backing o 

Set up matrix P^^ of Legendre polynomials for experimental 
angles • 

Do least squares calculation to find coefficients using the 
reciprocal of the estimated intensity variances from 11 as 
weights o 

Normalize coefficients and correct for the effect of the 

finite solid angle subtended by the scintillation crystal o 

These latter corrections are precomputed in an auxiliary 

41 

program using the expressions given by Rose« Estimated 

variances of the coefficient ratiosj A«/A^, etCo, are com^ 

^ ^ 61 

puted from the expression given by Reich et alo (After 

2 

computation variances were increased by a factor 1/^ in 

2 

all cases where the computed value of £ was less than unityo) 
Calculate anisotropyo 

2 

Compute the criterion £ <, This quantity, defined by Rose®s 

analysis, is a measure of non-statistical errors in the data. 

^ 2 

If such errors are eliminated the value of ^ serves to in“ 

dicate whether the appropriate number of pol)momials has been 

used in the correlation function expression, with values of 
2 

C greater than unity indicating an inappropriate choice. 

It is pertinent here that non^statistical errors have already 

been introduced into the problem in 11 above so that this 

criterion is valid only to the extent that an accurate esti“ 

mate of the true variance of the counts has been achieved. 

2 

The computed values of ^ in general are near unity and seem 
to verify that the estimated variances can be considered to 
be reasonably good. 
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go Selected values of the normalized calculated correlation 
function, and values of a comparison theoretical function 
if desired, are read out, and similarly normalized experi“ 
mental points are read outo (Experimental values are un~ 
corrected for solid angle, the correction being applied to 
the theoretical function for comparison with experimental 
results o ) 

ho Iterate for each gamma-ray energy o 

io Reduce number of polynomials if orders higher than have 

been used and iterate solutiono (The program also calculates 
coefficients for the cosine series: 



jo Return to main routine » 

16 o Repeat entire analysis for as many sets of data as indicated » 

17 o Combine coefficients to get weighted mean values and Chi square 
values if more than one set of data is analyzed o Read out these values o 
18. Combine intensities at each angle from all sets of data and re- 
analyze for coefficients using this combined data set. 

C. Glossary of names of selected variables and constants appearing 
in program SPECFIT. An asterisk indicates an input quantity. 

A(I,J) The matric product of the response matrix with its 




m 



) 



1 



transpose, weighted with the WT values. 



AINV(I, J) 



The inverse of A(I,J). 



ANISO(I) 



Anisotropy for Legendre series solution (1=1), 
and for cosine series solution (I =: 2) . 



ARGP(I) 



Vector containing cosines of the angles B3(I). 



*B(I) 



A vector representing the product of the inverse 
of the response matrix and the observed spectrum. 



*BACK(I) 



A number representing background counts to be 
subtracted from calibration curve I. 



BACKGND(I) 



The calculated background spectrum to be subtracted 
from an input spectrum. 



*BCHAN 



The lower of the channel pair described under TCHAN 
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FIN(I,J,K) 


Storage for FINTENS(I), data set Jj run K« 


FMU(I) 


Mass absorption coefficient in copper for EMU(I) 
Mev gamma ray Units in cm^/gm» 


*HCM 


Distance from source to crystal face in centimeters « 


HOLD (I) 


A temporary storage vector » 


LPKS 


Same as NPKS but not including any straight line 
component 0 (See NBAOPT) 


*MAXCHAN 


An upper channel limit in loading CALB datao Between 
MAXCHAN and MINT, CALB(I) values are filled in with 
the value of CALB (MAXCHAN) « 


*MESH 


The number of channels in input spectrum » 


*MINT 


Number of channels in calibration spectra « 


*MINTB 


Number of background channels « 


MPKS 


Initial value of NPKS for restoring the value. 


MPKS2 


Initial value of LPKS for restoring the value. 


*NBAOPT 


Background option selector. 

NBAOPT =0 No background read in nor subtracted 
from SPEC 

NBAOPT = 1 Background spectrum is read in and 

SPEC is corrected for this background 

NBAOPT = 2 Same as NBAOPT = 1 and in addition a 
straight line is included in the 
response matrix. 


*NCUT 


A channel number below which the graph output 
values are divided by a scale factor SCALE » 


NERFL 


An output of matrix inversion subroutine. A value 
of 1 indicates matrix is non singular, and 2 indicates 
a singular matrix. 


*NHIEST 


The upper channel of the range over which a fit is 
made in SPECFIT. 


*NLOWEST 


The lower channel of the range over which a fit is 
made in SPECFIT. 


*NOCALB 


Number of calibration spectra. 


*NPHOTO 


The photograph number of the observed spectrum; for 
identification only. 
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*NPKS 


The number of spectral components with which the 
observed spectrum is to fitted o Program will increase 
this value by one if option is chosen to include a 
straight line in R (See NBAOPT) » 


*NRUNS 


Number of runs in a set of data* 


NP 


Number of Legendre polynomials to be fitted in first 
pass through LEASTSQ* 


*NPRINT 


Print suppression index 

NPRINT = 0 All printouts and graphs are made 

NPRINT = 1 Suppress printouts of spectrum values and 
cosine series coefficients 

NPRINT = 2 Same as 1 and also suppress graph outputs 

NPRINT = 3 Same as 2 with briefer more compact arrange 
ment of correlation coefficients 

NPRINT = 4 Same as 3 but graphs plotted 


*NSETS 


Number of sets of data. 


OMEGA(I,J) 


Used for temporary storage in SPECFIT and diagonal 
elements used for weights in LEASTSQo 


*PKSEP(I) 


The number of channels corresponding to 1,022 Mev 
for energy calibration of CALB curves. Taken as 
channel distance from photo peak to pair peak for 
higher energy calibration gammas. 


*PHOT(I) 


The channel of calibration curve I in which the 
photopeak occurs, (May be an interpolation between 
two integral channels numbers). 


*PKEN(I) 


Energy in Mev of a spectral component. 


POLY(l.J) 


The matrix of Legendre polynomials of order 2(J-1) 
evaluated at angle B3(I)« 


R(I.J) 


The response matrix. 


RED 


A normalizing factor for 3** x 3** unshielded Nal(Tl) 
crystal representing the probability that one photon 
at the source will produce one count in the photopeak 
channe 1 . 


*RES 


The proton energy in kev at the resonance being 
investigated (for identification only). 


SAN (I) 


Estimate of standard deviation of ANISO(I) , 


SAVE (I, J) 


An auxiliary matrix. 
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SBB(I) 


Coefficient for renormalizing input spectrum with 
counts from an auxiliary monitor (optional; if not 
entered values are unity) « 


* SCALE 


A scale factor for graph output <, See NCUT. 


SCOEFT(I) 


The standard deviation in COEFT(I)<. 


SIGCTS(I) 


The calculated error in COUNTS(I)- 


SIGTENS(I) 


The error for calculated spectral intensity 
INTENS(I) 0 


SIN(I,J,K,) 


Storage for errors in FIN(IsJ,K) 


*SPEC(I) 


The vector representing the counts in the observed 
spectrum. 


STCCljJ^K) 


Storage for accumulating calculated correlation 
coefficients of order I, for gamma ray J and 
data set K, 


STORE(I,J) 


A temporary storage matrix for A while it is being 
inverted. 


STSC(I,J,K) 


Standard deviation of STC(I,J,K) , 


STW(I,J) 


Temporary storage. 


STSW(I,J) 


Temporary storage. 


SUM7(I) 


Temporary storage. 


*TBACK 


Duration of background run in units of seconds x 60 


TC(I) 


The weighted mean value of correlation coefficient 
of order I . 


*TCHAN 


The upper of a pair of input spectrum channels to be 
used for sensing energy calibration of the input spectrum. 


TEMP (I) 


Temporary storage. 


*TOPEN 


Energy of TCHAN in Mev, 


TS(I) 


Standard deviation of TC(I) 


*TX60SEC 


Duration of a run in units of seconds x 60, 


WT(I) 


The weight assigned to channel I in least squares 
solution for intensities. 


XI(I) 


The matric product of POLY and Bl weighted with the 
diagonal elements of OMEGA, 
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ivsr 




m - 

I 






XS(I) 



A vector representing the energy in Mev of each 
channel of the input spectrum, 

Z(I) A temporary storage vector, 

*ZERO The number which if any diagonal element of matrix 

A is smaller than, A will be considered singular. 
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Do DATA INPUT SEQUENCE FOR PROGRAM SPECFIT 



Input 


Number 
of sets 


Number of 
cards in set 


Input name 


Remarks 


1 


1 


1 


NSETS, NPRINT 


4 


2 


1 


1 


EL, ECC, HCM 




3 


1 


1 


NOCALB, MINT 




4 


NOCALB 


1 


BACK(J), PHOT(J), ENCALB(J), 
MAXCHAN, PKSEP(J) 








MAXCHAN/7 


CALB(I,J) 




5 


1 


1 


NPKS, NBAOPT 




6 


1 


NPKS 


PKEN(I) 




7 


1 


1 


BSL, BEN, BPHOT, MINTB, 
TBACK 


1 


8 


1 


MINTB/7 


BTGT(I) 




9 


1 


1 


NRUNS, CTSNORA, NLOWEST, 
NHIEST, NCUT, SCALE, ZERO 


2 


10 


1 


1 


RES, NPHOTO 


3 


11 


1 


1 


MESH, TCHAN, TOPEN, BCHAN, 
BOTEN, B3(NRUNDEX), T60XSEC 




12 


1 


1 


CTSNORB 




13 


1 


MESH/7 


SPEC (I) 




14 


1 


4 


Graph data cards 


5 



REMARKS o 

1.) Read in only for NBAOPT greater than lo 

2o) Input resumes here for each successive data set* 

3o) Input resumes here for each successive run within a data seto 

4o) All loading formats are multiples of Fl0o3 and IlO except that ZERO 
is ElOol and TRACK is F20.3<, 

5o) Blank cards may be loaded if graph option not usedo 
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E. Description of program CONTOUR 



As described in appendix 1 the correlation function coefficients for 
(p,y ) and (p,xY) intermediate unobserved reactions for spin 1/2 target 
nuclides can be described as a product of, at most, three functions of the 

f orm ; _ 2 

^ik ^ *^ik k ^ '^ik O k 

®ik "" 



1 +^; 



where i denotes the order of the coefficient and k is the index of 

the mixing parameter^ . For computation the substitution ^ = tan Y 

60 

is commonly used, and leads to g = OJ + p sin ( 2 Y + JT ) 

a + c 

with (X = — ^ — 








1/2 



y = tan ^ - - ) 

s = b/ ] b I 

where indices have been dropped for clarity. 

This program computes sets of possible values of two mixing para* 
meters for various fixed values of correlation coefficients. Using a 
graph plotter output medium available through a locally developed sub* 
routine, these sets of values are plotted to form a set of contours 
along which the value of the correlation coefficient is constant. 

Steps in the computation are as follows: 

lo Read in data for one possible spin sequence to be examined. 
Compute values of a, P and JT for the lowest order coefficient 
which occurs. If three mixing parameters occur simultaneously 3 
read out the expression for the fixed g as a function of the 
fixed mixing parameter. (All loading formats are standardized 
to be F10.3 and 110 for floating and fixed point variables 
respectively. ) 

2. Find and read out the maximum and minimum possible values of 
the coefficient and the corresponding values of the Y*s. 
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3, 


Fix the coefficient at an initial value within its range o 
Assign values to one of the Y*s in equal increments through- 
out the range -90° to +90° and solve for the corresponding 
values of the remaining Y, which in general will be double 
valued over the region selected o 


4, 


Read out pairs of Y values for the selected coefficient 
value both to the print medium and to the graph plotter. 


5. 


Repeat steps three and four for successively larger values 
of the coefficient until a specified number of contours has 
been plotted. 


6. 


Repeat steps two through five until contour plots for co- 
efficients of all orders permitted by the given spin sequence 
have been produced. 



F. Input variables for CONTOUR 



NSETS 


Number of spin sequences for which plots are to be 
made , 


JFIND 


Index of the mixing parameter to be solved for, (Y 
values on graph) , 

1 = proton partial wave mixing 

2 = unobserved gamma multipolarity mixing 

3 = observed gamma multipolarity mixing 

4 = channel spin mixing 


JSTEP 


Index of mixing parameter assigned successive values 
(x values on graph) See JFIND, 


JFIX 


Index of mixing parameter to be held constant when 
three parameters occur simultaneously, (Set to five 
for (p, ^ ) reaction where only two parameters occur,) 


IMAX 


Index specifying the highest order coefficient occurring 
in the spin sequence, (Order = 2# IMAX) 


JMAX 


Largest value which the index of mixing ratio assumes 
for use in reading input data. (Usually set at four,) 


NCONTS 


Number of contours to be plotted on each graph. 


NPTS 


Number of pairs of points to be computed for each branch 
of one contour. 


DELTA (JFIX) 


Fixed value assigned to the third mixing parameter during 
a calculation. 
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RHO Energy dependent Coulomb phase factor to be multiplied 

into the b coefficients for proton^par tial^wave 
mixing, 

A(I,J) The a, b, and c factors occurring in the expression for 

B(I,J) go I denotes the coefficient order (order - 2l) and 

C(I,J) J the mixing parameter index. See FIND, 
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